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The model of self-organizing Eulerian walkers is numerically investigated on the square lattice. The 
critical exponents for the distribution of a number of steps (77) and visited sites (t s ) characterizing 
the process of transformation from one recurrent configuration to another are calculated using the 
finite-size scaling analysis. Two different kinds of dynamical rules are considered. The results of 
simulations show that both the versions of the model belong to the same class of universality with 
the critical exponents n = t s — 1.75 ± 0.1. 
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To understand the nature and physical origins of the 
Self-Organized Criticality (SOC) [jj, a number of differ- 
ent models have been suggested in the last few years. 
Typically, these models evolve according to the pre- 
scribed dynamical rules into the SOC state where they 
show spatio-temporal self-similarity and are character- 
ized by long range correlations. In the critical state, the 
systems pass from one stable configuration to another 
through the avalanches which play an essential role in 
the organization of the dynamical criticality. 

In this letter, we numerically investigate critical dy- 
namics of a new cellular automaton model recently pro- 
posed by Priezzhev, Dhar et al || . This model has some 
common features with the well-known Abelian Sandpilc 
Model (ASM) |],[| but differs in rules which govern the 
motion of particles on the lattice. 

The model of the self-organizing Eulerian walkers on 
the square lattice is defined as follows. We associate with 
each site i of the two-dimensional L x L lattice an arrow 
directed up, right, down or left from i. We start with 
an arbitrary initial configuration of arrows on the lat- 
tice. Initially, we drop a particle on the randomly chosen 
site i. The succeeding steps the particle performs are 
determined by the following rules: 

i) the particle coming to a site j turns the arrow clock- 
wise by the right angle, 

ii) then makes a step along the new direction of the 
arrow to the neighbor site, 

iii) if the new direction points out the lattice, the par- 
ticle leaves the system. 

These rules are applied until the particle eventually 
leaves the lattice. Then, we go on by adding a new par- 
ticle and so on. In this type of dynamics, the movement 
of the particle affects the medium and in turn is affected 
by the medium. 

On the lattice with closed boundary conditions, the 
particle never leaves the system and finally gets into a 
limit cycle in which it passes each bond in both directions 
only once. The walks of this type are known as Euler 
circuits ||. 

Let us consider the lattice with open boundary condi- 
tions. The set of bonds marked by arrows form a graph 



G. Adding the particles followed by their movement 
through the lattice changes the configuration of arrows 
and organizes the system into the SOC state This 
critical state is a collection of recurrent configurations 
of a Markov process because the motion of a succeding 
particle depends only on the final configuration of ar- 
rows produced bu its predecessor. Each element of this 
recurrent set can be obtained from the previous one by 
adding a particle followed by its movement and leaving 
the lattice. It turns out that for any recurrent configura- 
tion, the graph G is a spanning tree on the lattice. Thus, 
the set of recurrent configurations is in the one-to-one 
correspondence with spanning trees ^ . 

At each intermediate step, the moving particle can de- 
stroy a spanning tree and form a loop of arrows. At this 
moment the system leaves the recurrent set. Eventually, 
after finite number of steps the particle reconstructs the 
structure of a spanning tree. The interval between the 
destruction and restoration of the spanning tree can be 
called an avalanche of cyclicity. During the walk of the 
particle, the system passes several times from one recur- 
rent configuration to another through the avalanches of 
cyclicity. This process is similar to the avalanche dy- 
namics of sand in the ASM where avalanches also re- 
construct the recurrent configurations (or spanning trees 
as the correspondence between the recurrent set of ASM 
and spanning trees does exist as well). 

At the beginning of an avalanche of cyclicity the last 
turned arrow closes a loop. Then, the trajectory of the 
particle covers the interior of the loop. During this walk 
each inner arrow turns four times, whereas the arrows 
forming the loop turn in such a way that the direction 
of the loop gets reversed. It might occur that a moving 
particle, after escaping from one closed loop, may form a 
new loop. Due to this structure of the walk, the number 
of steps in an avalanche is equal to /c(4n +1), where k 
is the number of loops constituting the avalanche and 
n = 0, 1, . . . It is possible to prove that the avalanche 
may consist of only one or two loops. This fact explains 
the line doubling in the distribution of the steps in the 
avalanche of cyclicity (Fig.l). 

To investigate the avalanche process in the models of 
self-organizing Eulerian walkers, we studied them nu- 
merically with high statistics. For each distribution of 
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avalanches we considered up to 30 • 10 6 events on the 
square lattices of linear size L from 120 to 400. Simula- 
tions always started from the regular initial configuration 
in which all arrows were directed up. 

In Fig. 2, we present the double logarithmic plot of the 
distribution P(s) of the number of visited sites in the 
avalanche for the lattice size L = 400. The analysis of 
the data shows that this distribution obeys the power law 



P(s) 



(1) 



To estimate the critical exponents, we have performed 
the finite-size scaling analysis |^J^] assuming the distri- 
bution functions scale with the lattice size L as 



P{x,L) = L-Pf{x/L») 



(2) 



where f(y) is a universal scaling function and and v 
are critical exponents which describe the scaling of the 
distribution function. 

To reduce the fluctuations of the data, we inte- 
grated each distribution over exponentially increasing bin 
lengths. For the integrated bin distribution we have § 



D(s) = J P{x) dx 



,-(t,-1) 



(3) 



Plotting D(s,L)L^ 3 versus sL~ Vb on a double loga- 
rithmic scale, as is shown in Fig. 3 for the different lattice 
sizes L, we obtained that the best data collapse corre- 
sponds to (3 8 = 1.5 ± 0.05, v s = 2.0 ± 0.05 (Fig.4). The 
scaling relation for the critical exponents t s — (3 s /v s -\- 1 
gives the value r s — 1.75 ± 0.05. 

In the same way, we investigated the distribution P(l) 
of the number of steps performed by the particle in the 
avalanche for the different lattice sizes L. There is an 
explicit power law behavior in these distributions (Fig.l) 



p(0 ~ r 



(4) 



with the finite size cutoff. We applied the finite-size scal- 
ing analysis to the integrated distributions and obtained 
77 = 1.7 ± 0.05 from the best data collapse. 

We also investigated a slightly modified model. The 
difference from the previous one is in the order of turns of 
the arrow. In the case when the turns form the sequence 
(up-down- left-right-up), we found a similar power law for 
avalanche distributions. 

To find critical exponents of this power law from finite- 
size scaling analysis, we integrated again these distribu- 
tions over exponentially increasing bin lengths. The crit- 
ical exponents 77 = 1.72 ± 0.05 and r s = 1.8 ± 0.05 have 
been obtained from the best data collapse for the distri- 
bution of steps and visited sites, respectively. 

In conclusion, we numerically investigated the model 
of self-organizing Eulerian walkers on the square lattice. 
The dynamics of the model organizes the medium of the 
system and builds up spatio-temporal complexity. We 
obtained explicit power law distributions in two slightly 



different versions of the model. We calculated the crit- 
ical exponents for the distribution of a number of vis- 
ited sites (r s ) and number of steps (77) in avalanches of 
cyclicity. These exponents are equal within a small un- 
certainty We argue that the critical exponents for these 
models within small errors belongs to the same class of 
universality and have a surprisingly large value 1.75±0.1 
in comparison with the known exponent for the ASM 
(r = 5/4) [§. 
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FIG.l The distribution P{1) of the number of steps in avalanches on the square lattice of the linear size L = 400. 
FIG. 2 The distribution P(s) of the number of visited sites in avalanches on the square lattice of the linear size L = 400. 
FIG. 3 The integrated distributions D(s) for the eight lattice sizes with L ranging from 120 to 400. 
FIG. 4 The finite-size scaling for the integrated distributions D(s). 
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FIG. 1. The distribution P(l) of the number of steps in avalanches on the square lattice of the linear size L = 400. 
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FIG. 2. The distribution P(s) of the number of visited sites in avalanches on the square lattice of the linear size L = 400. 




FIG. 3. The integrated distributions D(s) for the eight lattice sizes with L ranging from 120 to 400. 
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FIG. 4. The finite-size scaling for the integrated distributions D(s). 



